Dynamic scaling in the spatial distribution of persistent sites 
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The spatial distribution of persistent (unvisited) sites in one dimensional yl + A model 
is studied. The 'empty interval distribution' n{k,t), which is the probability that two consecutive 
persistent sites are separated by distance k at time t, is investigated in detail. It is found that at late 
times this distribution has the dynamical scaling form n{k, t) ^ t~^ k~^ f{k/t^). The new exponents 
r and z change with the initial particle density no, and are related to the persistence exponent 9 
through the scaling relation z{2 — r) — 6. We show by rigorous analytic arguments that for all no, 
1 < r < 2, which is confirmed by numerical results. 
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First passage problems in non-equilibrium systems un- 
dergoing time evolution has become an important field of 
research lately with the discovery of persistence. Persis- 
tence probability in general is defined as follows: Given 
a stochastic variable 4'{t) which fiuctuates about a mean 
value, say zero, what is the probability P{ti, that (j){t) 
does not change sign throughout the time interval [ii, i2]- 
For a large class of physical systems, persistence shows a 
power-law decay P{ti,t2) ~ (^2/^1)"^ for ^2^*1, with 
a non-trivial persistence exponent [p]-p2[ which is, in 
general, unrelated to other known static and dynamic 
exponents. 

Let us consider spatially extended systems with a 
stochastic field 0(x, t) at each lattice site x, the time 
evolution of which is coupled to that of its neighbouring 
sites. (t){x,t) could be, for instance, an Ising spin [0J^, a 
phase ordering field , a diffusing field or the height 
of a fluctuating interface In such cases, the sys- 

tem gets broken up into domains of persistent and non- 
persistent sites in course of time. In d = 1, this reduces 
to a set of disjoint persistent and non-persistent clusters 
appearing alternately. As persistence decays with time, 
the persistent clusters shrink in size and hence, their sep- 
aration grows. The following questions arise naturally in 
this context, which we address here: (i) How are the per- 
sistent clusters distributed in space at a given time? (ii) 
How does their average separation grow with time? 

In one dimension, the zeroes of the stochastic field can 
be viewed as a set of particles, moving about in the lat- 
tice, annihilating each other when two of them meet. 
When a particle moves across a lattice site for the first 
time, the field there flips sign, and the site becomes non- 
persistent. If each particle is assumed to perform purely 
diffusive motion, this reduces to the well-known reaction- 
diffusion model A-l-A —> 0, with appropriate initial condi- 
tions. The simplest case is random initial distribution of 
particles, with average density uq, for which P{t) ^ t^^ 
with 9 = 3/8 independent of no We investigate 
spatial ordering of persistent sites in this simple model. 

Our study is centered around the Empty Interval Dis- 



tribution n{k, t) — the probability that two randomly 
chosen consecutive persistent sites are separated by dis- 
tance k at time t. This distribution is analogous to the 
well-studied Inter-Particle Distribution Function (IPDF) 
in diffusion- reaction systems jl^ . Our numerical results 
show that n{k, t) has a non-trivial dynamic scaling form 
with power-law decay in k and t, characterized by expo- 
nents T and 9 respectively. The power law decay is valid 
for k <^ L{t), where L{t) ~ is a new dynamic length 
scale, which may be interpreted as the average separa- 
tion between persistent regions. The three exponents are 
connected by the scaling relation z{2 — t) = 9. Although 
the persistence exponent 9 is universal for this model, 
we flnd that r and z do change with the initial particle 
density. We give rigorous analytical arguments on the 
bounds of t, showing that 1 < t < 2 for all values of tiq. 
Power-law decay of n{k,t) in fc is a consequence of spa- 
tial correlations — a random distribution of sites would 
correspond to exponential decay. 

Our numerical simulation is done on a 1-d lattice of 
size N = 10'* with periodic boundary conditions. Par- 
ticles are initially distributed at random on the lattice, 
and their positions are sequentially updated — each par- 
ticle was made to move one step in either direction with 
probability 1/2. When two particles came on top of each 
other, both vanished instantaneously. The time evolu- 
tion is done up to 12000 Monte-Carlo steps (1 MC step is 
counted after all the particles in the lattice were touched 
once). All simulations are repeated for three different val- 
ues of initial density, no = 0.2, 0.5 and 0.8. The results 
are averaged over 500 different initial configurations. 

We observe that at large times t and fc 1, n(fc, t) ~ 
k~'^ for k <^ L{t). Here L{t) is a cut-off length scale 
that grows with time. In Fig. 1, we present the data for 
— 0.5 and for three values of time. The same data as 
presented in Fig. 2 shows that for each fc, nk{t) ^ t^'^ 
for all k <^ L{t). (It will be shown later that lu ~ 9). 
Similar power-law decay in fc and t has been observed 
for other values of uq also. These observations are fairly 
well represented by the following dynamic scaling form 



1 



for n{k,t), for late times and large enough k. 

n{k,t)^t-^k--f{k/L{t)) (1) 

where the scaling function f{x) ~ 1 for a; ^ 1 and 
decreases faster than any power of x for x 3> 1. 
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Fig. 1 Log-log plot of n{k, t) with k illustrating the power 
law decay for small fc, which crosses over to faster decay at 
large k. The data is presented for t = 2000, 4000 and 10000 
(top to bottom). The vertical separation between the curves 
has been enhanced for clarity. The initial density no is 0.5. 
All the three straight lines have slope r ~ 1.27. 
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Fig. 2 Log-log plot of n{k, t) with t for fe = 1 (top) and 
fc = 10 (bottom). Initial density no is 0.5. The asymptotic 
slope of the lines gives uj ~ 0.37. 

The exponents appearing in Eq. |^ are not all indepen- 
dent. The moments of the distribution are useful in de- 
riving the scaling relations between them. The m-th mo- 
ment is Im{t) = X]fe ^'"'^e^' ^ n{s,t)s"^ds. From 
the definition of n(fc, t), one can easily see that 

Io{t) = P{t) ^ t-'' ; h{t) = N (2) 

The average separation between persistent sites is 
given by l2{t)/h{t) - L[t). In Fig. (3) we have L{t) 



plotted against t on a logarithmic scale for three values 
of riQ. We find that L{t) diverges with time as 

L(t) ^ (3) 

where z is a new dynamic exponent. 
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Fig. 3 The average separation between persistent sites 
L{t) grows as a power of time t. The three logarithmic plots 
correspond to no = 0.8, 0.5 and 0.2 (top to bottom). The 
lines are visibly getting flatter with increasing no, indicating 
the decrease in the exponent z. 

The scaling relations between the exponents are ob- 
tained by making use of the conditions in Eq. ^ First 
of all, we show that only r < 2 is physically reasonable. 
For, if T > 2, Ii{t) ^ and from the second part of 
Eq. H it follows that uj — 0. But since uj > 9 for reasons 
of convergence, we get < which is absurd. So we con- 
clude that r < 2. In this case, Ii{t) ~ ^-'-'+^(2-^)^ which 
according to Eq. ^ imply that 

z(2 - r) = w (4) 

Another set of scaling relations can be derived using 
the condition on /o(t) in Eq. ^ Combined with Eq. |^, 
this gives 

z = e ] uj = e{2-T) if T < 1 (5) 

LU = e ; z{2-t) =e if r > 1 (6) 

We present a summary of our numerical results in Ta- 
ble 1. ft is easily seen that for all values of tiq, z > 9, 
T > 1 and the scaling relations in Eq. |^ are satisfied 
within numerical errors. Moreover, the exponents z and 
r show a consistent decrease with increasing no — they 
are non- universal, unlike 9. We now present an intuitive 
argument which accounts for these observations fairly rig- 
orously. 

As persistence decays, the non-persistent regions grow 
in time (the length scale of which is set by L{t) in Eq. 
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while the clusters of persistent sites shrink in size and 
eventually disappear. Let p{l,t) be the number of per- 
sistent clusters of size I, at time t. The total number 
of persistent sites at time t is J^i^Pi^^^) = ^^'^ 
the total number of such clusters is Nc{t) = J2iPi^T^)- 
The latter is related to n{k, t) through the exact relation 

Nc{t) = X]fcL2 ^(^' ^) ^ ^(^) ^ "'(-'^' '^^^ average size 
of a cluster at time t is 



m_ 

Nc{t) 



= 1- 



Pit) 



(7) 
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From Eq. 0, n{l,t) - t-'^ and since P{t) ~ t 
have = [l ~ jt^^^^^''] where 7 is a numerical con- 
stant. Since ui > 0, l{t) is a constant for late times. Now, 
if a; > 6*, l{t) = 1 strictly; only iiu — 9 any other value is 
possible. We argue for the latter case as follows. When 
clusters of persistent sites shrink in size, the depletion 
happens at the two ends of the cluster, independent of 
its size. Let the average decrease in the size of a clus- 
ter over time t be £,{t). Clusters of initial size I > ^{t) 
shrink to size / — ^(i) after time t, while those with length 
I C(^) disappear. It follows that 



(8) 



Here, p{l, 0) = nQ^l — no)' since the initial distribution 
of particles is done at random with probability no. Af- 
ter substitution in Eq. ^ we find that the time evolution 
of the cluster size distribution has the extremely simple 
form p{l, t) = e-^«(*V(', 0) where A = -ln(l - no). This 
result is also supported by simulations (Fig. 4). Con- 
sequently, the average cluster size J{t) = 1(0) = I/uq. 
This implies uj — 6 from our arguments following Eq. 
and thus validates Eq. ^. Furthermore, since P{t) ~ 
we find S,{t) ~ jln t at large t. It follows that a per- 
sistent cluster of initial size L has an average life-time 
tl ~ exp{^L) for large L. The exponential dependence 
of the life-time of the cluster on its size reflects the slow 
algebraic decay of persistence. 

Our argument can be extended to show why r and 
z are possibly non-universal. First of all, the exponent 
relation lu = 9 makes it possible to write the formal re- 
lation n{l,t) = g{T,no)P{t) Combined with Eq. 
and using the result l{t) — 1/no, we obtain the relation 
g{T,no) + no = 1, which expresses implicitly the depen- 
dence of T on riQ. For instance, if Eq. ^ were exact for 
all values of fc, then P{t) « n(s, t)ds = "J^'*-* so that 
g(r, no) ~ T — 1, from which it follows that r ~ 2 — no. 
This result, although not exact, is consistent with the 
bounds 1 < r < 2, and appears to be valid in the high 
density limit no — > 1, as indicated by the numerical val- 
ues in Table I. 
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Fig. 4 The cluster size distribution p{l, t) decays expo- 
nentially with I at all times. The figure shows p(l,t) plotted 
against I for t = 20, 40 and 120 (top to bottom) and no — 0.2. 
The vertical separations have been enhanced for clarity. 

In summary, we have shown that the spatial distri- 
bution of persistent clusters in one dimension exhibits 
rich dynamic scaling characterised by two new exponents. 
We have given rigorous arguments on the bounds and 
universality properties of these exponents, which is well- 
supported by numerics. Interestingly, the normalized size 
distribution of persistent clusters was found to be inde- 
pendent of time. 

Our work is the first study that brings out the non- 
trivial features in the spatial distribution of persistent 
sites in a one dimensional model. The dynamic scaling 
form in Eq. |^ is by no means specific to the model stud- 
ied, and we have observed similar forms in other one di- 
mensional systems — diffusion equation and kinetic Ising 
model, for example. Similar scaling in size distribution 
has been observed in entirely different contexts also — 
for instance, diffusion-limited cluster aggregation [T5|-[T7|] 
and diffusion- limited deposition |l8) . The feature that is 
common to all these processes is the irreversible coales- 
cence of clusters (empty intervals in our case) . 

We are grateful to Satya Majumdar for discussions and 
for pointing out the similarities to aggregation models. 
G. M thanks B. Derrida and P. R thanks D. Stauffer for 
critical reading of the manuscript and valuable sugges- 
tions. 



no 


0.2 


0.5 


0.8 


e 


0.3718 ± O.QQQl 


0.3769 ± 0.0000 


0.3729 ± 0.0001 




0.3763 ± 0.0001 


0.3767 ± 0.0000 


0.3710 ± 0.0002 


z 


0.5766 ± 0.0005 


0.5107 ± 0.0003 


0.4392 ± 0.0005 


r 


1.3502 ± 0.0262 


1.2596 ± 0.0357 


1.1277 ± 0.0383 



TABLE I. Exponents 6, lo, z and r as measured from the 
simulations, for three values of initial density no. The exact 
value of is 0.375, independent of no. The uj values presented 
correspond to = 1 and the r values tot = 12000. The errors 
given are only statistical. 
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